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Abstract 

The correlator of singlet axial-vector and vector currents in the external elec- 
tromagnetic field is studied within the instanton liquid model of QCD vacuum. In 
the chiral limit we calculate the longitudinal and transversal Wj! with respect 
to axial- vector index invariant amplitudes at arbitrary momentum transfer q. It is 
' demonstrated how the anomalous longitudinal part of the correlator is renormalized 

Q^ ■ at low momenta due to the presence of the V a (1) anomaly. 

in 

\n ■ 1 Introduction 
o 

Consideration of the axial-vector A and vector V current-current correlator in the soft 
external electromagnetic field V is an important part of the calculations of the compli- 



Ph. cated light-by-light scattering amplitude related to the problem of accurate computation 
of higher order hadronic contributions to muon anomalous magnetic moment 1 . In this spe- 
cific kinematics when one photon (V) with momentum q<i = q is virtual and another one 
(V) with momentum q\ represents the external electromagnetic field and can be regarded 
as a real photon with the vanishingly small momentum qi only two invariant functions 
survive in linear in small q\ approximation. It is convenient to parameterize the VAV 
correlator (Fig. ^) in terms of longitudinal Wl and transversal Wt with respect to axial 
current index Lorentz invariant amplitudes 



T^x(qi, 92) = [ w t {q 2 ) pi** - <?2<?i<?2^a + q\q{q2 e P^v) - w L (q 2 ) qlq{qle mav \ . 

(1) 

Both Lorentz structures are transversal with respect to vector current, q^T v \ = 0. As for 
the axial current, the first structure is transversal with respect to q\ while the second one 
is longitudinal and thus anomalous. The appearance of the longitudinal structure is the 
consequence of the Adler-Bell-Jackiw axial anomaly 

For the nonsinglet axial current A^ there are no perturbative [5] and nonperturba- 
tive [H] corrections to the axial anomaly and, as consequence, the invariant function w L 

1 See, e.g., PG] and references therein. 
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Figure 1: Diagrammatic representation of the triangle diagram in the local perturbative 
theory (a); and in the instanton model with dressed quark lines and full quark-current 
vertices (b). 

remains intact when interaction with gluons is taken into account. It was shown in [7j 

(3) 

that in the nonsinglet channel the transversal structure Wj, is also free from perturba- 
tive corrections. Nonperturbative nonrenormalization of the nonsinglet longitudinal part 
follows from the 't Hooft consistency condition jH], i.e. the exact quark-hadron duality 
realized as a correspondence between the infrared singularity of the quark triangle and 
the massless pion pole in terms of hadrons. However, for the singlet axial current due 
to the gluonic Ua (1) anomaly there is no massless state even in the chiral limit. Instead, 
the massive 77' meson appears. So, one expects nonperturbative renormalization of the 
singlet anomalous amplitude at momenta below rf mass. 

In the previous work [S] we analyzed in the framework of the instanton liquid model 
jHj the nonperturbative properties of the nonsinglet triangle diagram in the kinematics 
specified above. We demonstrated how the anomalous structure w L is saturated within 
the instanton liquid model. We also calculated the transversal invariant function w T 
at arbitrary q and show that within the instanton model at large q 2 there are no power 
corrections to this structure. The nonperturbative corrections to at large q 2 have 
exponentially decreasing behavior related to the short distance properties of the instanton 
nonlocality in the QCD vacuum. 

The present work is devoted to the study of the Ua (1) anomaly effect on the singlet 
VAV correlator within the instanton liquid model. We calculate in the chiral limit the 
longitudinal and transversal invariant functions at arbitrary momentum transfer 
q and demonstrate how the singlet anomalous part of the correlator is renormal- 
ized at low momenta due to presence of the Ua (1) anomaly. 

2 The structure of VAV correlator in perturbative 
approach 

The amplitude for the triangle diagram can be written as a correlator of the axial current 
jl and two vector currents j u and j M (Fig. QJ) 

T^x = - J d 4 *dV^-^ (0| T{j v (x)j„(y)jl(0)}\0) , (2) 
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where for light u and d quarks one has in the local theory 



Jm = q It^Vq, fx = Q IxlsAq , 

the quark field q^ has color (2) and flavor (/) indices, A(°> = I, = r 3 are flavor matrix 

of the axial current, and V = V = \ {\ + T z) are the charge matrices, with the tilted 
current being for the soft momentum photon vertex. 

In the local perturbative theory the one-loop result (Fig. for the invariant functions 
wt and wl for space-like momenta q (q 2 > 0) is 



w ^ 1o °p = 2^- 1oop 



2N c Tr(AVv) f 2 , (3) 



where N c is the color number, and for light quark masses one takes rrif = m u « m^. In 
the chiral limit, rrif = 0, one gets the result 

w L (q 2 ) = 2w T (q 2 ) = ^Tr (AVv) . (4) 



3 The instanton effective quark model 

To study nonperturbative effects in the triangle amplitude T^ u \ at low and intermediate 
momenta one can use the framework of the effective approach based on the representation 
of QCD vacuum as an ensemble of strong vacuum fluctuations of gluon field, instantons. 
Spontaneous breaking the chiral symmetry and dynamical generation of a momentum- 
dependent quark mass are naturally explained within the instanton liquid model. The 
instanton fluctuations characterize nonlocal properties of the QCD vacuum [H H ITT | IT2]. 
The interaction of light u, d quarks in the instanton vacuum can be described in terms of 
the effective 't Hooft four-quark action with nonlocal kernel induced by quark zero modes 
in the instanton field. The gauged version of the model [T3*l IT4"1 ITSj meets the symmetry 
properties with respect to external gauge fields, and corresponding vertices satisfy the 
Ward-Takahashi identities. 

In the framework of this effective model the nonsinglet V and A current-current cor- 
relators, the vector Adler function, the pion transition form factor have been calculated 
for arbitrary current virtualities in ^HIGE1E|- In the same model the topological suscep- 
tibility of the QCD vacuum which is reduced to the singlet A current-current correlator 
has been considered in jTHl HE] • 

The spin-flavor structure of the nonlocal chirally invariant interaction of soft quarks 
is given by the matrix products 2 

G (1 ® 1 + i^r a <g> i 75 r a ) , G' (r a g> r a + i 7s <g> i lb ) , (5) 

where G and G' are 4-quark couplings in iso-triplet and singlet channels. For the inter- 
action in the form of 't Hooft determinant one has the relation G' = —G. In general due 
to repulsion in the singlet channel the relation G' < G is required. 

Within the gauged instanton model the dressed quark propagator in the chiral limit, 
S(p), is defined as 

£r\p) = iv- M(p 2 ), (6) 

2 The explicit calculations below are performed in SUf{2) sector of the model. 
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with the momentum-dependent quark mass 



M(p 2 ) = M q f(p 2 ) (7) 
found as the solution of the gap equation 



d A k , 2n2 ,M(k 



M(p 2 ) = AG P N f N c f(p 2 ) J ^-J\k 



2\ 



(2vr) 4J v D{k 2 ) ' 
where we denote 

D(k 2 ) = k 2 + M 2 (k). 

The constant M q = M(0) in (J7J) is determined dynamically from Eq. (JHJ) and the function 
f(p) defines the nonlocal kernel of the four-quark interaction. Within the instanton model 
f(p) describing the momentum distribution of quarks in the nonperturbative QCD vacuum 
is expressed through the quark zero mode function. It is implied in jSJ HU H2] that 
the quark zero mode in the instanton field is taken in the axial gauge when the gauge 
dependent dynamical quark mass is defined. In particular it means that f(p) for large 
arguments decreases like some exponential in p 2 . To make numerics simpler we shell use 
the Gaussian form 

f(p) = exp ( -p 2 /A 2 ) , (9) 

where the parameter A characterizes the size of nonlocal fluctuations in the QCD vacuum 
and it is proportional to the inverse average size of an instanton. 

The vector vertex following from the instanton model is (Fig- Et) 

r^k, k') = ^ + (k + k')^^, k'), (io) 

where M^(k, k') is the finite-difference derivative of the dynamical quark mass, q is the 
momentum corresponding to the current, and k (k') is the incoming (outgoing) momentum 
of the quark, k! = k + q. The finite-difference derivative of an arbitrary function F is 
defined as 

= 3£^. (id 

The nonlocal part of the vertex (jlOjl necessarily appears in order to fit the vector 
Ward-Takahashi identity. 

Within the chiral quark model based on the non-local structure of instanton 
vacuum jTTj the full singlet axial- vector vertex including local and nonlocal pieces is given 
by [IHl 

If»{k, q,k' = k + q)= 7m75 + 75 (fc + k%M q U ^I'Jjl^ 2 ~ ( 12 ) 

i <2V r n '\ f a \ G' 1 _ GJ PP (q 2 ) 
+ 7s^ 2M q f(k)f{k) 



where 



Y qJ v /J v ' Gl-G'J PP (q 2 ) 



_ ( 2 .. 8N C f d 4 k MM (/,-. •// • ) 

JMq ~ ~m% J (^j 1 r^rZ • (13) 
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Figure 2: Diagrammatic representation of the bare (a) and full (b) quark-current vertices. 



Here and below we use the notations 

k-q, 



k + 
M± 



k ; k - 
M(kl), 



kl + M 2 



±1 



k + k- 
f± 



(k+q) {k-q) 



f(k 



The vertex ()12|) takes into account the quark- ant iquark rescattering in the singlet 
axial channel (Fig. |2b)- For completeness we also present the vertex corresponding to the 
conserved iso-triplet axial- vector current which in the chiral limit is given by ^3] 



rf\k,k') = r 3 



lfll5 + 2 l5 q ^MJ(k)f(k') 



(14) 



+l 5 {k + fc%M< 



{f{k')-f{k)f 
k' 2 - k 2 



The isotriplet axial-vector vertex has a massless pole, q 2 = 0, that follows from the 
spontaneous breaking of the chiral symmetry in the limit of massless u and d quarks. 
Evidently, this pole corresponds to the massless Goldstone pion. 

The singlet current (fT2"|) does not contain massless pole due to presence of the Ua (1) 
anomaly. Indeed, as q 2 — > there is compensation between denominator and numerator 
in (H2J) 



1 - GJp P (q 2 
-q 2 



G 



M 2 q 



as q 



0. 



(15) 



where is the pion weak decay constant. In cancellation of the massless pole the gap 
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equation is used. Instead, the singlet current develops a pole at the 77'— meson mass 3 

1 -G'J PP (q 2 = -m 2 v ,) = 0, (16) 

thus solving the Ua{^) problem. Let us also remind that in the instanton chiral quark 
model the connection between the gluon and effective quark degrees of freedom is fixed 
by the gap equation. 

4 Singlet VAV correlator 

In the effective instanton-like model the nondiagonal correlator of vector current and 
singlet axial-vector current in the external electromagnetic field (VAV) is given by (Fig. 

d A k 

•Tr [r„ (k + q 1} k) S (k + 9l ) rf } (k + q u k - q 2 ) S (k - q 2 ) T u (k, k - q 2 ) S (k) 

where the quark propagator, the vector and the axial- vector vertices are defined by ©, 
([TUjl and JUJ), respectively. The structure of the vector vertex (fTU|) guarantees that the 
amplitude is transversal with respect to vector indices 

Tp,\(qi, g 2 )9i = f^ uX (qt, q 2 )q 2 = 

and the Lorentz structure of the amplitude is given by (JTJ). 

In jH] we found for the nonsinglet axial current the expressions for the longitudinal 
amplitude 

(3), 2< 2N, 1 



T^ x {q u q 2 ) = -2N c Tt (AW) / 7—4- (17) 



3 g 2 ' v ' 

and for the combination of invariant functions which shows up the nonperturbative dy- 
namics 

,2' 



w 



§(g 2 )^4 3) (? 2 )-24 3) (9 2 




v/M+M (1) (A; + ,A;_) - 2 (kq) M' + Vm \k+, k„) \, (19) 



3 L 

where prime means a derivative with respect to k 2 : M'(k 2 ) = dM(k 2 )/dk 2 . The result 
(J 18)) which is independent of the details of the nonlocal effective model is in agreement with 
the statement about absence of nonperturbative corrections to the nonsinglet longitudinal 
invariant function that follows from the 't Hooft duality arguments. 

The calculations of the singlet VAV correlator results in the following modification of 
the nonsinglet amplitudes 

wP(q 2 ) = 5 -W^(q 2 )+Aw^(q 2 ), (20) 
4V) = jjtl® (q 2 ) + Aw^ (q 2 ) , (21) 



3 See previous footnote. Also we neglect the effect of the axial-pseudoscalar mixing with the longitu- 
dinal component of the flavor singlet f± meson. 
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where 



Aw® (g 2 ) 



5N C 1 - G'/G f d 4 k ^M+M 



9g 2 1 - G'Jpp (g 2 ) 

'4 {kqf 
3 g 2 



7T 4 DlD- 



Mj 



- -Ml A; 2 
3 + x 



-M {1 \k + ,k_ 



(kg) 



(22) 



In jH] it was shown that in the chiral limit the nonsinglet transversal amplitude gets 
only exponentially suppressed at large momenta corrections. The reason is that the 
asymptotics of the amplitude (jl9j) is proportional to the vacuum nonlocality function 
/ (g) that is necessarily has exponentially decreasing asymptotics. The singlet amplitudes 
differ from the nonsinglet ones by the term (J22|) that also has exponentially suppressed 
large q 2 asymptotics. Thus, within the instanton model the singlet longitudinal and 
transversal parts have only exponentially suppressed at large q 2 corrections and all allowed 
by operator product expansion power corrections are canceled each other. 
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Figure 3: Normalized wl invariant 
function in the singlet case (solid line) 
and nonsinglet case (dashed line). 



Figure 4: Normalized wlt invariant 
function versus Q predicted by the in- 
stanton model in the singlet case (solid 
line) and isotriplet case (dashed line). 



Fig. 3 illustrates how the singlet longitudinal amplitude Hi® is renormalized at low 
momenta by the presence of the Ua (1) anomaly. The behavior of w^(q 2 ) is presented in 
Fig. 4. In both figures the corresponding results for the nonsinglet case are also shown. 
The values of the model parameters used in calculations are fixed earlier in jHJ EI] as 



M q = 0.24 GeV, 



A P = 1.11 GeV, G P = 27 A GeV" 



(23) 



The coupling G' is fixed by fitting the meson spectrum. Approximately one has G' ~ 0.1 G 
|13j . We also find numerical values of the invariant amplitudes at zero virtuality 



w f\ q 2 = 0) = 4.4 GeV~ 



wfltf = 0) = 0.6 GeV" 



(24) 



Precise form and even sign of wf^{q 2 ) strongly depend on the ratio of couplings G'/G 
and has to be defined in the calculations with more realistic choice of model parameters. 
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5 Conclusions 



In the framework of the instanton liquid model we have calculated for arbitrary momenta 
transfer the nondiagonal correlator of the singlet axial-vector and vector currents in the 
background of a soft vector field. For this specific kinematics we find that in the chiral 
limit the large momenta power corrections are absent for both longitudinal wl as well 
transversal wt invariant amplitudes. These amplitudes have very similar behavior and 
are corrected only by exponentially small terms which reflect the nonlocal structure of 
QCD vacuum. 

Within the instanton model the renormalization of the singlet longitudinal am- 
plitude occurring at low momenta due to the Ua (1) anomaly is demonstrated explicitly. 
In the nonsinglet case the behavior of wl and wt at low momenta is very different due to 
the contribution of the massless pion state. At the same time in the singlet case there is 
no massless state and the deflection of wl from 2wt amplitudes is rather small. 

The author is grateful to A. P. Bakulev, N. I. Kochelev, P. Kroll, S. V. Mikhailov, A. A. 
Pivovarov, O. V. Teryaev for helpful discussions on the subject of the present work. The 
author also thanks for partial support from the Russian Foundation for Basic Research 
projects nos. 03-02-17291, 04-02-16445 and the Heisenberg-Landau program. 

References 



[i] 

[2] 
[3] 
W 
[5] 
[6] 

[7 
[8 
[9 
[10 

[11 

[12 
[13 



S. Groote, J. G. Kdrner, A. A. Pivovarov, Eur. Phys. J. C 24, 393 (2002). 

A. Czarnecki, W. J. Marciano and A. Vainshtein, Phys. Rev. D 67 (2003) 073006. 

S. L. Adler, Phys. Rev. 177, 2426 (1969). 

J. S. Bell and R. Jackiw, Nuovo Cim. A 60, 47 (1969). 

S. L. Adler and W. A. Bardeen, Phys. Rev. 182, 1517 (1969). 

G. 't Hooft, in Recent Developments In Gauge Theories, Eds. G. 't Hooft et at, 
(Plenum Press, New York, 1980). 

A. Vainshtein, Phys. Lett. B 569 (2003) 187. 

A.E. Dorokhov, arXiv:hep-ph/0505007| 



See for review, e.g., T. Schafer, E.V. Shuryak, Rev. Mod. Phys. 70 (1998) 323. 

S.V. Mikhailov, A.V. Radyushkin, Sov. J. Nucl. Phys. 49 (1989) 494 [Yad. Fiz. 49 
(1988) 794]; S. V. Mikhailov and A. V. Radyushkin, Phys. Rev. D 45 (1992) 1754. 

A. E. Dorokhov, S. V. Esaibegian and S. V. Mikhailov, Phys. Rev. D 56 (1997) 4062; 
A.E. Dorokhov, S.V. Esaibegyan, A. E. Maximov, S.V. Mikhailov, Eur. Phys. J. C 
13 (2000) 331. 

A. E. Dorokhov and L. Tomio, Phys. Rev. D 62 (2000) 014016. 
R.S. Plant and M. C. Birse, Nucl. Phys. A 628 (1998) 607. 



S 



[14] I. V. Anikin, A. E. Dorokhov and L. Tomio, Phys. Part. Nucl. 31 (2000) 509 [Fiz. 
Elem. Chast. Atom. Yadra 31 (2000) 1023]. 

[15] A. E. Dorokhov, W. Broniowski, Eur. Phys. J. C 32 (2003) 79. 

[16] A. E. Dorokhov, Phys. Rev. D 70 (2004) 094011. 

[17] A. E. Dorokhov, JETP Letters, 77 (2003) 63 [Pisma ZHETF, 77 (2003) 68]; 
I. V. Anikin, A. E. Dorokhov and L. Tomio, Phys. Lett. B 475 (2000) 361. 

[18] A. E. Dorokhov, Phys. Part. Nucl. Lett. 1 (2004) 240. 



9 



